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I. INTRODUCTION 

A mixing angle 6*13 in the lepton sector has not been measured yet, and we have only an 
upper bound obtained by the CHOOZ reactor experiment [1]. Recently reactor experiments 
have attracted much attention again as a possibility to measure ^13 [2-11]. One of the 
advantages of the reactor measurement of ^^13 is possible resolution of the 6*23 ambiguity [3, 
12, 13]. For the reactor measurement to resolve the 623 ambiguity, however, the error in 
the reactor measurement of ^13 has to be relatively small. The error in the measurement of 
sin^ 2^13 is almost independent of the central value of sin^ 2^13 [3] , so it is basically the same 
as the sensitivity to sin^2^i3. To resolve the ^23 ambiguity by the reactor measurement, 
therefore, the sensitivity to sin^ 2^13 has to be reasonably small. Thus estimation of the 
sensitivity to sin^ 26*13 in the reactor measurement is important, and for the purpose of 
designing reactor experiments in the future where the statistical error becomes negligibly 
small due to large volume of the detectors, it is useful to know what gives "the systematic 
limit", i.e. the sensitivity to sin^ 2^^i3 in the limit of infinite statistics. 

Refs. [2-4] discussed the idea of the near-far detector complex which improves the sensi- 
tivity to sin^ 2^13, and Ref. [14] discussed analytically what gives the dominant contribution 
to the systematic limit on sin^ 2^13 by a rate analysis. In this paper we discuss the system- 
atic limit using an analytical approach which was described in detail in [14] . The advantage 
of the analytical treatment is that we can easily see which systematic error gives dominant 
contribution to the systematic limit on sin^ 26*13 in the rate as well as spectrum analysis. 

Throughout this paper we discuss the case with a single reactor and two detectors for 
simplicity in the limit of infinite statistics. Also for simplicity, we assume in this paper that 
the near and far detectors are identical and have the same sizes of systematic errors. In 
Sect. 2 we review the results of a rate analysis to this system and show that the dominant 
contribution to the sensitivity to sin^ 26*13 comes from the uncorrelated error between 
the detector. We also derive lower bound on the sensitivity, which is basically proportional 
to o"u. In Sect. 3 we perform a spectrum analysis on the same reactor system and show 
the dependence of the contribution from each systematic error on the sensitivity. In the 
appendices we give some details on how to derive the analytic results used in the main text. 

II. SYSTEMATIC LIMIT OF sin^ 26*13 BY A RATE ANALYSIS AND ITS LOWER 
BOUND [14, 15] 

Let us first review the sensitivity by a rate analysis in the limit of infinite statistics. Let 
and be the number of events measured at the near and far detectors, and be the 
theoretical predictions. We must consider the correlation of errors between the detectors. Uc 
and (7u indicate the correlated and uncorrelated systematic errors in the number of events, 
respectively. Then is given by [14] 
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where e{E), Y(E), (t{E) stand for the detection efficiency, the neutrino yield at the reactor, 
and the cross section, respectively. The numerical value of D{L) is plotted in Fig. 1 as a 
function of L, where the reference value for lAmJgl is 2.2xlO~^eV^ which was obtained by a 
combined analysis [16] of the data of the atmospheric neutrino and K2K experiments. Here 
we assume (Ju = C(0.1)% and (Tc = C(l)%. (4(7^ + 2al)~^ is much smaller than (2cr^)~^ 
while D{L-f) + D{L]<f) and D{L-p) — D{L^) are generally comparable in magnitude, so the 
first term in Eq. (2) can be ignored.^ Hence is given approximately by 

^^^^r.',eJ^M^p^. (3) 

The hypothesis of no oscillation is excluded at the 90%CL if is larger than 2.7, which 

corresponds to the value at the 90%CL for one degree of freedom. This implies that the 
systematic limit on sin^ 26*13 at the 90%CL, or the sensitivity in the limit of infinite statistics, 
is given by 

, sys only 



V limit ^(^p) _ D(Ln) ^ > 

To optimize (sin^ 2^13) j'j^""^^, therefore, we have to minimize D{L^) = (sin^ (Aml^L^/AE)) 
and maximize D[L-p) = (sin^ (Aml^Lp / 4E)) . Since the possible maximum value of D{Lp) — 
L'(Ln) is 0.82, which is obtained with Lp = 2.0km and Lp = for lAm^g] = 2.2 x 10"^^ 
the lower bound of (sin^ 2^13)^^°'^'^ in this case is estimated as: 

lower bound of (sin^ 29,^) ZT^ ^^^^ = 2.8 a^. (5) 

Eq. (5) indicates that the sensitivity is at best sin^ 2^13 ~ 0.016 if we adopt the reference 
value (Tu = 0.6% assumed in Ref. [3]. To reach this hmit of the sensitivity, we have to 

set Lp as close to 2.0km as possible, set Ln as close to as possible. For the case with 
finite statistical errors, Lp which is slightly smaller than 2km is appropriate because it gives 
smaller statistical errors. 



III. SYSTEMATIC LIMIT OF sin^ 26*13 BY A SPECTRAL ANALYSIS 

In the previous section we have seen that the sensitivity to sin^ 2^13 is governed by the 
uncorrelated systematic error of the detectors in the rate analysis. In this section we will 
examine what happens if we use spectral information also. 



^ In the region of very large | Am^gl, the first term dominates the sensitivity because of D{Lp) ~ D{Ln) — 
0.5. 
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We introduce the systematic errors in almost the same way as in [4]. ^ is given by 
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Here, mf is the number of events to be measured at the near [A = N) and far [A = F) 
for the i-th energy bin with the neutrino oscillation, and tf is the theoretical prediction 
without the oscillation, {afy is the uncorrelated error which consists of the statistical plus 
uncorrelated bin-to-bin systematic error: 

where a^^ is the uncorrelated bin-to-bin systematic error. For simplicity we assume that 
sizes of the bin-to-bin uncorrelated systematic errors of the detectors and of the flux are 
independent of the energy. This may not be the case in practical situations, but our analytic 
discussions will be illuminative to see how each systematic error affects the sensitivity, a is 
a variable which corresponds to a common overall normalization error ctdb for the number 
of events. «^ {A = N, F) is a variable which introduces the detector-specific uncertainties 
(TdB of the near and far detectors, (i = 1, • • • , n) is a variable for an uncertainty a^h of 
the theoretical prediction for each energy bin which is uncorrelated between different energy 
bins.^ a^i {A = N. F) is a variable which introduces an energy calibration uncertainty a^ai 
and comes in the theoretical prediction in the form of (1 + a^J-E instead of the observed 
energy E. Thus, the deviation vf (divided by the expected number of events) from the 
theoretical prediction tf due to this uncertainty can be written as 



lim 
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where Np is the number of target protons in the detector, T denotes the exposure time, and 
La is the basehne for the detector A. We have used the definition 
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(8) 



Our strategy is to assume no oscillafAon in the theoretical prediction tf, to plug the number 
of events with oscillation in mf, and to sec the sensitivity by looking at the value of x^. 
Since is quadratical in the variables a, a^, ct^, ct^i, we can minimize with respect to 



^ Here we follow the notation for the systematic errors in Ref. [17]. The first suffix of a stands for the 
property for the systematic error with respect to the detectors while the second is with respect to bins, 
and capital (small) letter stands for a correlated (uncorrelated) systematic error. The correspondence for 
the notation in Ref. [4] is as follows: (T„ = crab) fft = crdB, Cshape = CDbi <^a = cdb- 
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these variables in Eq. (6) exactly.^ Hereafter, we take the hmit of infinite statistics, i.e., 
1 / tf ^ in order to see how systematic errors affect the sensitivity to sin^2^i3; We refer 
to the sensitivity in this hmit as "the systematic hmit". Since we assume no osciUation for 
the theoretical prediction tf, the quantity vf in Eq. (7) is independent of A{= N, F): 

(9) 

After some calculations (see Appendix A for details), (6) is rewritten as 
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where we have introduced the variables 



u-'-iD^ + D^) , (10) 
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and are the orthonormal eigenvectors of some nxn unitary matrices (see Appendix A) . 
a' ^ is defined by h^~^ = a"^^ and (Jcai while A^^^ is defined by b^~^^ = cr^g + a^B and (Tcai; 
a' ^ (A^^') vanishes for b^^^ = (6*^+^ = 0) or deal = 0, and then the second (fifth) term in 
Eq. (10) degenerate to the first (fourth) term. The first three terms in Eq. (10) correspond 
to the comparisons of the numbers of events at the two detectors because these terms does 
not include + and the correlated errors between the detectors {a-ah and (Jdb)- We 



^ In principle we could take a different convention, such as multiplying the measured numbers mf by the 

uncertainty (1 + a + + a^), etc., and in fact it is done in some references. With such a convention, 
it becomes complicated to work with an analytical approach because is not a Gaussian with respect 
to mf^ , which includes oscillation parameters, after the minimizations of a's. Note that = 2.7 gives a 
90%CL bound for a Gaussian. However, the difference between such a convention and ours affects only 
the higher orders in cr^'s. The conclusion on the sensitivity with such a convention should coincide with 
ours numerically. 
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see that the first and fourth terms in Eq. (10) correspond to the "pure" spectral analysis 
because these terms are free from CTcai and correlated errors among bins (udB and Udb); 
Other four terms include the effect of the rate analysis. 

Here we will take the following assumptions for the systematic errors: 
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(12) 

Roughly speaking, the errors must satisfy the relations ct^ ~ a^g + (T%f,/n and a1 ~ (j^g + 
a^^/n, where cTu and cTc are the errors appear in the rate analysis (2). Among the reactor 
experiments in the past, Bugey [18] seems to be the only experiment in which the energy 
spectrum analysis was performed with the identical detectors. We have learned [19] that a^h 
in the Bugey experiment was of order of 0.5%. We have checked numerically that each factor 

n—2 

- E W ■ ± D"")] \ - [Ci • {D^ ± D"")] \ - [ut' ■ {D^ ± ^^)] ' is approximately 

independent of the number n of bins for n ^ 8,'^ and the maximum value of these factors 
range from 0.4 to 2 (See Table 1). In most of the analyses in the present paper we take the 
number n of bins n = 16 and the energy interval 2.8MeV< <7.8MeV. The optimized 
baselines Lp and L^v do depend on the energy interval of the analysis because different 
energy interval gives different average of the neutrino energy. We adopt the energy interval 
2.8McV< Ely <7.8MeV because (i) the same energy interval was used in the analysis of 
the Bugey experiment, and (ii) the value of A_ /n becomes approximately independent of 
the number n of bins, (ii) is not necessarily the case, for example, if we take the energy 
interval 1.8MeV< Ej, <7.8MeV because in this case the value of \v\ becomes so large near 
the threshold energy E^, ~ 1.8MeV that /n would depend on n. In fact, we can neglect 
such an energy region near the threshold because the region has only small number of events 
and does not affect the value of so much. 

Prom the assumed values in Eq. (12) it follows that the last three terms in Eq. (10) are neg- 
hgible because the coefficients in front of these three factors are much smaller than those in 

^ n-2 

front of the first three and because numerical values of the factors — 
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1 



n I 



^(±) 



are all comparable, as can be seen from 

Table I. Furthermore, it turns out that the third term in Eq. (10) is negligible compared to 
the second term. This is because 



^ We found from numerical calculations that the results with different numbers (n = 2™, to = 3, • • • ,6) of 
bins do not change much. Fig. 2 gives the dependence of the first three quantities in Eq. (10) on the 
number n of bins iov Lp = L and Ln = 0, and it shows that the behaviors for n > 8 are almost the same. 
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are satisfied for our reference values (See Appendix C) and because the values of 
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This is the main result of this paper and is to be compared with the result (3) by the rate anal- 
ysis. Once the baselines Lp and L^r are given, then (13) contains only (t\^/ y/n and cr^g as the 
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parameters because the coefficients — ^ 
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are almost independent of n. These coefficients seem to give good measures for the potential 
of an experiment in spectral analysis almost independently of n and the sizes of errors; The 
larger values of these measures means the better setup in principle. 

From numerical calculations we find that the dependence of the first term in Eq. (13) 
on Lp and is quite different from that of the second term and the values of 

1 2 
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Ln in Figs. 4 and 5 
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of the first term in Eq. (13) has the maxi- 
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of the second 



is extremely small for these baselines.^ Theoretically, therefore, the first term is optimized 
for Lp =10. 6km and LAr=8.4km. In this case the sensitivity to sin^ 2^13 in the limit of 
infinite statistics is given by 



(sm 2^i3)ii^it 



X^ 90%CL ^ 2 ^db 



0.37 

for Lp — 10.6km and Ln — 8.4km. 
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If CTdh/ is smaller than Uu, the sensitivity (14) of the spectral analysis is better than (5) 
of the rate analysis. It should be noted, however, that this sensitivity is attained only if the 
statistical error in each bin becomes negligible compared to the systematic errors which are 
already assumed to be smaller than C(0.1)%. In order for this to happen at Lp =10km, 
we would need at least 100 kton-yr even at the Kashiwazaki-Kariwa nuclear power plant, so 
the optimization of the first term in (13), namely the optimization of the spectral analysis, 
is not reahstic for the investigation of the oscillation of |Am^3| = O{10~^)eV^. 



The values of the two optimized baselines for (1/n) 



-,(-) 



{D — D ) depend on the number n of bins, 



and the present values are obtained for n = 16. This is because the more number of bins we have, the 
more local maxima and minima wc could observe, and in order to see more local maxima and minima we 
should have longer baselines. Even if we take n > 16, however, this set of the baselines Lp =10. 6km and 
Ljv=8.4km gives the local maximum and can be regarded approximately as the optimized set. 
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Next, wc consider to optimize the the second term in Eq. (13), i.e. the term which roughly 
corresponds to the rate analysis. This is achieved at Lp =2. 1km and Ljv=Okm. For these 
values of the baselines, Eq. (13) becomes numerically 

2 sin^ 2^13 / 3.7 X 10-3 4.1 x 10-^ \ 

for Lp — 2.1km and Ljv — 0km. 

If o^y^/n » (l/100)(T^g then the second term on the right-hand side in Eq. (15) is dominant. 
In this case, the sensitivity is not so much different from that in rate analysis because the 
second term is understood as a rate-like one. On the other hand, if o"db/^ {^I^^^)^ab 
then the first term in Eq. (15) dominates the sensitivity, and it is given by 



'4 



Since the coefficient 4.1 x lO--*^ of the second term in (15) is much larger than the coefficient 
3.7 X 10-3 of the first term, (Xdb/v^ must be much smaller than in order to use (16) and 
to improve the sensitivity significantly beyond the one in the rate analysis . If the setup 
is optimized for the rate analysis, therefore, it is not easy for the spectral analysis to give 
us much better sensitivity than what the rate analysis does. On the other hand, since the 
difference between the coefficients of the first and second terms in (15) can be smaller for 
the setup which deviates from the optimal one for the rate analysis, the spectral analysis 
could compensate the loss of the sensitivity due to the deviation. 

To illustrate the usefulness of the formula (13), let us see if it reproduces the result 
in [4], where the number n of the energy bins was assumed to be 62, the energy interval was 
1.8MeV< Eu <8.0MeV, Li7=1.7km, Ljv=0.17km, and the reference values are |Am|]^| = 
3.0 X lO^^eV^, crdB=0.6%, and crdb=0.1% (the most optimistic case in [4]). We found from 
numerical calculations that if we perform the spectrum analysis with the energy interval 
1.8MeV< Ey <8.0MeV, the number n=62 of bins, the baselines Lp^l.Tkm and L;v=0.17km, 
then Eq. (13) is modified for |Am|i| = 3.0 x lO-^eV^ as 



2 ^ sin^2g,3 / 6.8X10-2 3.3 x 10-^ 
^ - 2 \ ay 62 + ai/62 + 0.95a^B ' ^ ^ 

for Lp — 1.7km and Ljv = 0.17km, 



where we have used the fact v^/n= 1.2 x 10^ and v- h/n = -7.7 for 1.8MeV< E^ <8.0MeV 
and n—62 to derive the factor 0.95. Note that the ratio of the coefficient 0.068 of the first 
term to that 0.33 of the second term is bigger than that in Eq. (15) because LF=1.7km is 
longer than the optimized baseline L ~1.4km for the rate analysis in the case of lAmg^^l = 

3.0 X lO-^eV^. This can be seen from Fig. 2, where the value of - ^ ^uf^ ■ {D^ - D^)^ 
increases for L larger than the optimized baseline (~2km) for the rate analysis, i.e. the 



baseline which optimizes — 

n 



1 r r ^ ^ ^ 1 2 



in the case of |Am|i| = 2.2 x 10 ^eV^. 
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Substituting the reference values (7(16=0.6%, (Tdb=0.1% in Eq. (17) and the value 2.7 of 
at the 90%CL for one degree of freedom, we have 

2.7 ^ (sin-2g,3);^r^ / 6.8 x 10'^ x 62 ^ 0.33 



2 V 10-6 0.95 X (6 X 10-3)2 

, (sin2 2^^i3)r-r'' 

~ 4.2 X 10^ X ^ 

2 

which gives 



(sin2 2^i3)r.r'' ^ \/— X 10"' ~ 1 X 10-3. 

V ^"^/ limit V 4 2 

This is almost consistent with the result in [4]. From this calculation it is obvious that such 
a good sensitivity was obtained in [4] because the uncorrelated bin-to-bin systematic error 
(7db was assumed to be very small and the number of bins was large. 



IV. OSCILLATION EXPERIMENTS FOR Am^ ~ C»(1)EV2 

In the previous section, in the case of reactor experiments to measure sin^ 2^^i3 for Am^ ~ 
(9(10^'^)eV^, wc have seen that the optimal baseline lengths for the spectral analysis are very 
long and it is not a realistic idea to accumulate large number of events at such baselines. If 
we want to observe neutrino oscillations for Am^ ~ (9(l)eV^, however, the optimal setup 
for the spectral analysis becomes realistic. With the standard three flavor scenario, there 
is little motivation to look for neutrino oscillations for Am^ ~ (9(l)eV^, but if we consider 
more general frameworks with sterile neutrinos then oscillations with Am^ ~ 0(l)eV^ are 
still possible. In this section we briefly discuss such a possibility for the sake of completeness. 

In ref. [20] it was shown that the scheme with two sterile neutrinos is still consistent 
with all the data, and the best fit values for l\m\^ and sin^ 2^i4 are 0.9eV^ and 0.058. 
Am^=0.9cV^ is the value for which the hypothesis for neutrino oscillations is excluded only 
weakly by the Bugey experiment [18]. If the (3-1-2) -scheme is realized by Nature, therefore, 
we should be able to observe neutrino oscillation at Am^ ~0.9eV^. By rescaling the basehnes 
used in (14) by 2.2 x 10"3eV^/0.9eV^ ~ 1/410, we see in the limit of infinite statistics that 
the optimum set of baselines for Am^ = 0.9eV^ is 

Lp — 26m, Ljv = 21m. 

If we put the two detectors with this set of basehnes, then in the limit of infinite statistics, 
assuming the number of bins n — 25, the bin-to-bin uncorrelated systematic error (7db=0.5% 
and the mixing angle sin^ — 0.058, we could have 

, sin^2e 0.37 n , , 

~ ~ 6 X 10^ fl8) 

~ 2 (5 X 10-3)2/25 - ^ ■ ^-'^^ 

Eq. (18) would imply that neutrino oscillation for Am^ ~0.9eV^ can be estabhshed with 
tremendous significance. This result is to be compared to the rate analysis with (3), where 
the set of optimum basehnes for Am^ ~0.9eV^ would be 

Lp = 5m, = Om. (19) 
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Assuming that these basehnes are technically possible, at sin^ 1Q — 0.058 would be 

2 sin^2e 0.82 , 

X 7 ^ - 4 X 10\ (20) 

^ 2 (0.6 X 10-2)2 ^ > 

From practical point of view, however, Eq. (19) would imply that both the detectors be 
placed inside of the reactor and Eq. (19) is not realistic. The value of would be smaller 
than Eq. (20) and much smaller than Eq. (18). Thus, the spectral analysis is much more 
advantageous than the rate analysis in search of neutrino oscillations for l\.vr? ~ C(l)eV2. 



V. DISCUSSION AND CONCLUSION 

We applied our analytical method to estimate the sensitivity to sin2 Id in the limit of 
infinite statistics, which we refer to as "the systematic limit", in the spectral analysis at 
a reactor neutrino oscillation experiment with one reactor and two detectors. In a rate 
analysis, the systematic limit (4) is set by the the uncorrelated systematic error of the 
detectors, and is given by 2.8(Tu [14]. In a spectrum analysis, we derived an analogous formula 
for ^2 which gives the systematic limit. The formula for contains two terms which could 
potentially dominate (13). One is the contribution from the bin-to-bin correlated detector 
specific systematic error (JdB and the other is the one from a^h / \Jn which is the uncorrelated 
bin-to-bin systematic error divided by the square-root of the number n of bins. The optimal 
setup for the spectral analysis is found to be unrealistic for t^w? — 0{W~^)eV^ because the 
optimal baselines are too long to accumulate large number of events. The optimization can 
be realistic, however, for Am^ = 0{l)eV^ to search for sterile neutrino oscillations. If the 
baselines are optimized for the rate analysis, the improvement of the sensitivity with the 
spectral analysis is not significant in most cases. On the other hand, the spectral analysis 
can be effective for the experiment whose setup deviates from the optimal one for the rate 
analysis. For example, the spectral analysis could compensate the disadvantage of the rate 
analysis in the KASKA experiment which has near detectors with relatively long basehnes 
(~ 300-400m), or that in the DCHOOZ experiment that has a far detector at a distance 
(^ 1km) shorter than the optimal one for the rate analysis. In any case, it is very important 
in the future reactor experiments to estimate a^B and especially a^h- 



APPENDIX A: DERIVATION OF THE COVARIANCE MATRIX 

After taking the hmit 1/tf ^ of infinite statistics, Eq. (6) becomes 
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and let us introduce a vector notation 
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Following the discussions in the Appendix A in [14], the matrix element of the covariance 



matrix can be obtained as the expectation value of yfyf: 
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where the normalization J\f is defined in such a way that (1) = 1. Prom a straightforward 
calculation, we have 

\yi yj / ^ (7db + (7db + () + Oij (Job + d Vj a^^■^. 

Thus we have 

7N \ / jjN 



= ( r ^ r nP~'{ If ) - sin^ 2^,3 {D^^,D^n P~' 



D 



F 



where the covariance matrix p is defined by 
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Here /„ is an n x n unit matrix, and G„ are n x n matrices defined by 

1 ••• 1 
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/ vf V1V2 
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where Vj (j = I,-- - ,n) is defined by Eqs. (7), and (9). Note that the covariance matrix 
does not include oscillation parameters but only errors. Diagonalization of the covariance 
matrix is useful to see which errors dominate x^. 
Now (Al) can be cast into a block diagonal by 



M N 
N M 



III -I 



2 V In In 



M + N 
M-N 



In In 
'In In 
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We obtain 

o sin^ 2^13 
X = — ^ — 



M + N 
M -N 



-1 



We can formally diagonalize the symmetric matrices M + N and M — N assuming the 
existence of {n — 2) orthonormal eigenvectors which are orthogonal to v and /i = (1, • • • , 1) . 
This is described in Appendix B. Now M+N and M — N are diagonalized by the orthogonal 



matrices C/^^' and U 



r(-) 



M + N = U'^'V'^'U'^' ^, 
M -N = U^'^V'-'^U^'^ ^, 



(A6) 



where 



P^*^ = diag {\f\--- ,A^*') 



and X^^^ are the eigenvalues of M ± A?". Prom (A5) and (A6) we obtain 



X 



sin^ 2^13 



X 



P'^^ 



V 



(-) 



sin^2^,3 J^l<'-(^^ + ^'^) 



2 



A 



(+) 



A 



(-) 



where u^'' are the orthonormal eigenvectors which appear in U^"^^ : 



Putting 



we have from (A2) 



U = yu^ , • • • , w„ j . 



M + N = a^^In + b'-^^Hn + cGn, 



(+) 



(-) 

a = cr; 



MB) 



c = cr: 



cal' 
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We will see in Appendix B that A'"*"' are given by 



(+) 

a 






= 1,-- 


■ ,n 


-2) 


(+) 

a 




u 


— n — 


1) 




(+) 

a 


+Ar 


u 


— n) 




(-) 

a 


+ A" 


u 


= 1,- 


• ,n 


-2) 




u 


— n — 


1) 






+ A^-^ 


u 









where A^^^ are given by 



Hence we obtain the (10)- Appendix C we will see that A_,_ /n are approximately 
independent of the number n of bins for the energy interval 2.8MeV< E^, <7.8MeV. 



APPENDIX B: DIAGONALIZATION OF THE nxn MATRIX 



In this appendix we show how we diagonalize the matrix 

ain + bHn + cGn, 

where is the nxn unit matrix, and Hn, Gn and Vj are given by Eqs. (A3), (A4) and Eq. (7), 
respectively. Since a/„ only shifts each eigenvalue by a, we will discuss the diagonalization 
of hHn + cGn here. Let us introduce the notation 

..(■). 

Vi/ 

— * 

When h ■ V 0, which is satisfied in the present case for Eq. (7), we can always find in 
n-dimensional Euclidean space n — 2 orthonormal eigenvectors Uj {j — 1, - ■ • , n — 2) which 

— * 

are orthogonal to h and v, i.e., 

h-Uj^v-Uj^O (j = 1, • • • , n - 2). (Bl) 

Prom Eq. (Bl) we have 

HnUj = Uj ■ h 
GnUj = Uj ■ V 
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Eq. (B2) indicates that Uj {j — 1, - ■ ■ , n — 2) are the n — 2 eigenvectors of bHn + cGn with 
the eigenvalue 0: 

(bHr, + cGn)uj^O = ,n-2). 

It turns out that the remaining 2 eigenvectors can be constructed out of h and v. Straight- 
forward calculations show the following: 

— * — * 

Hnh — nh, 

— * — * 

Gnh — [v ■ h)v, 

— * — * 

HnV — {v • h)h, 
GnV = {if^)v. 



Hence we obtain 

(bHn + cGn) 



h \ _ f bn cv ■ h \ { h 
V j \bv-h cv^ 



This suggests that multiplication of the vectors ^/bh and ^/cv hj bH^ + cG„ give the 
symmetric matrix 

The symmetric matrix on the right-hand side can be diagonalized as 

hn Vbcv-h \ _ f cosijj — sin / ^ \ f cos^' sin'^ 
Vbcv-h cv^ )~\smip cosip ) \ A_ i I - sin V' cos V' 



where ijj is given by 



Vbcv -h ,^ 

tan2t^ = -- — — B4 

{bn-ctP)/2 ^ ^ 



and the eigenvalues A± are 



A± = ±W +(Vbcvh). (B5) 



Therefore, if we define 



wi A ^ / cosip simp \ f y/bh 
W2 J ~ \ —smilj cos J \ y/cv J ' 



then from Eq. (B3) we have 
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It is easy to show that the normaUzations of Wi and W2 are \/A+ and \/A_, respectively. 
We see that the {n — l)-th and the n-th orthonormal eigenvectors are given by 



Un-l 
Un. 



( 



cos t/j + y/cvsim/j 



(B6) 



with the eigenvalues A_ and A+. 

We note in passing that (B6) has the correct behaviors in the limit 6 -C c or 6 ^ c. In 
the limit 6 ^ c we have 



Un-l 



( h-{v-h)v/v^ \ 
n — {v ■ hy /iP 

V 



with 



while in the limit 6 ^ c we obtain 



with 



A_ 



n — 



cv 



J 



( V — {v ■ h)h/n ^ 



iP — (v ■ Ky/n 
h 



v 



[v-hf 



n 



bn 



\ 
J 



Putting everything together, we obtain 

C/^ {hHr, + cG„) U = diag(0, • • • , 0, A_, A+) 

with 



(B7) 



(B8) 



— * 

where lii, • • • ,1*71-2 are n — 2 orthonormal vectors which are orthogonal to v and /i, Un-i 
and Un are given by Eq. (B6), A_ and A+ are given by Eq. (B5). Adding the unit matrix 
ain to Eq. (B8), we finally obtain 

{aIn + hHn + cGn) U — diag(a, • • • , a, a + A_, a + A+). 
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APPENDIX C: THE PROPERTIES OF v 



From numerical calculations, we find for n ^ 8 and for 2.8MeV< <7.8MeV that 

1 



n 



• /i ~ —3.6 = independent of n, 

—v'^ ~ 24 = independent of n. 

n 



If we introduce the rescaled variables 



1 



n 



^ 1 - 
h = — ^/i, 

then t'^ and v ■ h are almost independent of the number n of bins. Furthermore, we rescale 
other quantities as follows: 



A 



± 



nA±, 



where A± is given by Eq. (B5). A± is approximately independent of n because v and v ■ h 
are: 



iT2 



A 



± 



b + cv 
2 

6 + 24c 



± 



± 



h — cv"^ 



he 



6 -24c 



+ 3.626c. 



(CI) 



Eq. (CI) suggests that we can expand A± in either 6 or c depending on whether 6 <^ c or 
6 3> c. In the present case, our reference values for the (— ) sector are 

c = = om%r). 

1 ^ ^ 2 1 2 

Fig. 3 shows the values of - ^ ^uf^ ■ {D^ ± D^) , - u^l^ ■ {D^ ± D^) and 



C ■ {D- ± D 



at Lp — L, Ljv = with different values of c/b for a^h — 0.6%, 

and we observe that the case of 6 ~ c is well approximated by the limit b <^ c. Therefore, 
from Eq. (B7) we have 



A_ 
A_ 



y2 



^2 
CV 



b 1 



24 



24c 



0.466 
24c 



With the notations used in the text, we have 



aLV^ - c^dB l-{v-hf/nv' ~0.46(j^B, 
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optimized Li?/km 


optimized Ljv/km 


optimized value 






2 


10.6 


8.4 


0.37 


1 

n 


^a^i . {D^ - D^) 


2 




2.1 





0.41 


n 




2 




3.5 





0.51 


n-2 


^) 


2 


11.5 


11.5 


0.44 


1 

n 




2 




2.1 


2.1 


1.6 


1 
n 




2 




3.5 


3.5 


2.1 



TABLE I: The values of each factor in the Eq. (10) in the spectrum analysis. The number of 
bins is taken to be n = 16 here, but the values of these factors are almost independent of n for 
n ^8. The energy interval adopted here is 2.8MeV< <7.8MeV, as in the case of the Bugey 
experiment. [18] 
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L/km 

FIG. 1: D{L) = (sin^ ) ^ function of L in the rate analysis. D{L) has its maximum value 

0.82 at L=2.0km, and approaches to its asymptotic value 0.5 as L oo. jAmfgl = 2.2 x 10~^eV^ 
is used as the reference value. 
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0.5 hcT§al»^u 



0.4 
0.3 
0.2 
0.1 






2000 4000 6000 

L/m 



8000 10000 



FIG. 2: The values of - ^ (^n" • D^^ (the dash-dotted hnes), - (^u~li ■ D^^ (the thick soUd 

. ,1 „ , O r, 



Unes) and — I n„ • -D 
n 



(the dashed hnes) for different numbers n of bins in the hmit a^^^ ^ a^, 
where Lp = L, Ln = 0. The latter two are almost saturated for n > 4, while the first quantity is 
saturated only for n > 16. |Amf3| is 2.2 x 10~^eV^. 
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(Un-1 -0)2/1 

(Un-1 -0)2/1 
(Un-1 -0)3/1 
("n-1-Di'/n 




2000 



(c/b=10"? 
(c/b=10-! 
(c/b=10° 
(c/b=10l 
(c/b=10^ 
(c/b=10"J 
(c/b=10-! 
(c/b=10° 
(c/b=10l 
(c/b=10^ 



4000 6000 
L/m 



8000 



10000 



n 

l(-) 



FIG. 3: The values of — {u„\ • D 



(the sohd lines) and — [u„ ^ ■ D'' 
n ^ 



(the dashed lines) for 



different values ofc/6 at Lp = L, Ln = 0, where b 



reference values in the main text give c/b^ ' ^ 1 and the two quantities in the figure for c/b^ ' ^ 1 



= (0.6%)2 and c = a^^i- The 

(-) 



are approximated by the limit c/ft' ' ^ 1. The number n of bins is 16 and lAmfgl = 2.2 x 10 ^eV^. 
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Ln /m 



n-2 

FIG. 4: The values of - ^ [-uj"' • {D^ - D^] 



value 0.37 at Li?=10.6km and LAr=8.4km. The number n of bins is 16 and |Am^ 



as a function of Lp and Ljv . It has the maximum 

2.2xl0^3eV2. 



131 
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FIG. 5: The values of 



n L 



D 



as a function of Lp and L 



value 0.41 at Lp=2.1km and L 



N'- 



=Okm. The number n of bins is 16 and jAmfg 



AT. It has the maximum 
2.2 X lO-^eV^. 
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